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We consider the physics and collider phenomenology of quirks that transform nontrivially under 
QCD color, SU(2)w as well as an SU(N) ic infracolor group. Our main motivation is to show that 
the recent Wjj excess observed by CDF naturally arises in quirky models. The basic pattern is that 
several different quirky states can be produced, some of which /3-decay during or after spin-down, 
leaving the lightest electrically neutral quirks to hadronize into a meson that subsequently decays into 
gluon jets. We analyze LEP II, Tevatron, UA2, and electroweak precision constraints, identifying 
the simplest viable models: scalar quirks ("squirks") transforming as color triplets, SU(2)w triplets 
and singlets, all with vanishing hypercharge. We calculate production cross sections, weak decay, 
spin-down, meson decay rates, and estimate efficiencies. The novel features of our quirky model 
includes: quirkonium decay proceeds into a pair of gluon jets, without a fe-jet component; there is 
essentially no associated Zjj or jjj signal; and there are potentially new (parameter-dependent) 
contributions to dijet production, multi-W production plus jets, W'y, 77 resonance signals, and 
monojet signals. There may be either underlying event from low energy QCD deposition resulting 
from quirky spin-down, and/or qualitatively modified event kinematics from infraglueball emission. 



I. INTRODUCTION 

CDF has recently reported a 4.1c excess in the Wjj 
event sample in 7.3 fb _1 of data for dijet invariant masses 
between 120 — 160 GeV p] [2]. The excess comprises of 
hundreds of events in the ljj+$T channel, corresponding 
to a sizeable cross section, a ~ few pb to account for the 
smaller leptonic branching fraction and the efficiency of 
detecting the signal. Several explanations have already 
appeared [31 0] , as well as detailed discussion of the size 
and shape of the Standard Model (SM) background [5]. 

In this paper we present an explanation of this ex- 
cess invoking pair production of "quirks" [7J , defined be- 
low, which subsequently undergo radiative energy loss 
and weak decay, finally annihilating into dijets. There 
are two main pathways that begin with quirk produc- 
tion and end with a Wjj signal consistent with the CDF 
analysis after their cuts: The first path consists of elec- 
troweak production of a heavy-light quirk pair, where the 
heavy quirk /3-decays into a light quirk, emitting a W , 
and then the remaining quirk-anti-quirk system settles 
into a quirkonium ground state. The quirkonium decays 
into gluon jets that reconstruct an invariant mass peak. 
The second path consists of strong production of a heavy 
quirk-anti-quirk pair, where both quirks /3-decay to their 
lighter electroweak partners. Even though two W bosons 
are emitted in this process, one W decay can be missed 
by the CDF analysis, particularly when the decay prod- 
ucts are light quark jets that have energies below the 
CDF jet energy cuts. Schematic diagrams of the produc- 
tion and decay of the heavy-light quirk system and the 
heavy-heavy quirk system is shown in Fig. [I] 

In this paper we consider "squirks" - the scalar ana- 
logues of "quirks" [51 [7J, which are new fields trans- 
forming under part of the SM gauge group along with a 
new strongly-coupled "infracolor" group SU(N) ic . The 
strong coupling scale of the infracolor group, A ic , is as- 
sumed to be much smaller than the masses of all quirky 




FIG. 1. Two classes of squirk pair production and /3-decay 
that can lead to a Wjj final state with a dijet invariant 
mass peak are shown in schematic diagrams. The top, 
"one-armed lobster" diagram corresponds to weak produc- 
tion pp — > W* — ¥ w±Wg + uqvJ^, following by charged squirk 
/3-decay u± — > W ± uq, and finally annihilation of the neutral 
quirks into a pair of gluons or infragluons (</). The bottom, 
"two-armed lobster" diagram corresponds to strong produc- 
tion of a heavy quirk-anti-quirk pair, following by /3-decay 
of both heavy squirks, and finally annihilation of the neutral 
quirks into gluons. The latter production can pass the CDF 
analysis with a modest efficiency when the decay products of 
one of the W's are not detected. Quirky spin-down radiation 
(into color gluons, infracolor glueballs, etc.) is not shown for 
either diagram. 



(or squirky) fields transforming under infracolor. Con- 
sequently, infracolor strings do not break, and quirk (or 
squirk) pairs remain in a bound state even when they 
are produced with high kinetic energies. This leads to 
several interesting collider physics, model building and 



dark matter applications [Trl21|. (Other work on hidden 
valley models can be found in [22l424j .) Certain kinds of 
quirks have already been searched for at the Tevatron by 
the DO collaboration [55]. Typically, squirks and quirks 
transform as a vector-like representation of the SM gauge 
group as well as a vector-like representation of infracolor, 
and thus can acquire any mass without any additional 
Higgs sector^ 

Here is a summary of our basic strategy to explain the 
Wjj signal, while avoiding the many constraints from 
LEP II, Tevatron, and UA2: 

1. We consider squirks (scalars), rather than quirks 
(fermions), mainly to employ a renormalizable op- 
erator to lead to mass splittings, while also min- 
imizing the electroweak precision contribution re- 
sulting from this splitting. 

2. The Wjj signal arises from a combination of the 
two distinct squirk pair production processes shown 
in Fig. [I] The weak decay, u± — » W ± u 0l could be 
2-body or 3-body, depending on the mass splitting 
between u± and uq. 

3. All of the squirks are color triplets under QCD. 
Hence, the leading visible annihilation channel of 
a squirk-anti-squirk bound state system is into gg. 
This provides our source of the jet-jet resonance, 
and consequently, we predict the flavor of the jets 
in the excess observed by CDF to be pure glue with 
no heavy flavor. 

4. The electrically-neutral squirks are mixtures among 
the T3 = component of an electroweak triplet as 
well as an electroweak singlet. Hence, they do not 
couple to the 7 or the Z. This implies the tree-level 
production cross section of the lightest squirks at 
LEP II vanishes. 

5. The electrically-charged squirks have masses larger 
than 100 GeV, and so are beyond the LEP II sen- 
sitivity. 

6. Light squirk pair production uqUq results in a dijet 
invariant mass around 150 — 160 GeV. The pro- 
duction cross section is far smaller than Tevatron's 
dijet sensitivity |27j . making it safe from Tevatron 
searches. The production cross section of this reso- 
nance at the CERN SppS {y/s = 630 GeV) collider 
is less than 1 pb per infracolor. The UA2 bound is 
roughly 90 pb [25], and thus does not restrict N ic . 

At this point we should emphasize that some aspects 
of squirk production and decay can be calculated or sim- 
ulated with standard collider tools, but some cannot and 
one must resort to estimates or parameterization of the 



1 Noteable exceptions, where quirks acquire masses through elec- 
troweak symmetry breaking, can be found in Refs. 16 20 . 
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various possibilities. What can be calculated unambigu- 
ously is i) squirk pair production (at leading order, far 
from the cross section threshold) , ii) the weak decay rate 
of squirks (from which we use to infer the weak decay 
of the squirky mesons), and iii) squirky meson annihila- 
tion rates into Standard Model fields. What cannot be 
calculated reliably is the dynamics of the "spin-down" 
process, as the high energy squirks shed their momen- 
tum to settle into a (near) ground state squirky meson, 
and the spectrum and content of the resulting radiation. 
An attempt at modeling this radiation is underway [21] . 
As a consequence, we do not attempt a full simulation 
of squirky production. Instead, we calculate quantities 
that are under control, with estimates of signal efficien- 
cies, and point out where simulations could improve our 
understanding. To deal with the remaining uncertainties 
due to non-perturbative infracolor dynamics we discuss 
some of the possible scenarios for this spin down and their 
effects on phenomenology in Sec. [IV| 

II. MODEL 

We extend the Standard Model to include one set of 
squirks, U and S with quantum numbers given in Ta- 
ble [I] The gauge-invariant operators involving these fields 
include the mass terms 

^M^U+hl^S (1) 

and quartic interactions including 

A 4 (S' t S*)(C/ t C/) + X V {H ] H){U ] U) + X S (H^H)(S^S) 
+ X Ui (U^U) 2 + A S4 (S t 5) 2 (2) 

as well as 

/s(.ff t 7- a .ff)(5 t I7 )+h.c.. (3) 

The operators proportional to Xu,s lead to shifts in the 
masses of U and S, as well as Higgs boson trilinear and 
quartic interactions with WlJ, S. The operators pro- 
portional to Xua,S4 will turn out to affect certain anni- 
hilation channels. Since our analysis does not depend on 
the existence of these interactions, we do not consider 
them further. 

The triplet field U can be written in terms of its isospin 
components 
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where we can rewrite the fields in the mass eigenstate 
basis with definite electric charges, 

u± = 4= (Ui T iU 2 ) (5) 

This is completely analogous to the rewriting of the usual 
triplet of SU(2) L of gauge bosons W 1 , 2 ,d, into W ± '°. 

After electroweak symmetry breaking, the Higgs dou- 
blet can be written as H = (v + h) exp[ill]/v2 in terms 
of the vacuum expectation value v = 246 GeV, the phys- 
ical Higgs boson h, and the nonlinear representation of 
the Goldstone bosons, contained within II. Expanding 
the operator Eq. ^ to second order in II, one can verify 
that only (v + h) enters the dynamics, and thus we can 
ignore the Goldstone interactions. 

The central implication of Eq. ([3| is that it causes 
the singlet to mix with the neutral component of the 
triplet U3, enlarging the mass splitting between them. 
The two mass eigenstates, which we call uq and U\ for 
the light and heavy states, have masses of 
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where i5 = v/ nv 2 /2. These neutral squirks are mixtures 
of the singlet and triplet 
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where we have used the notation Co = cos 9, sg = sin#, 
and the neutral state mixing angle is 

2<S 2 

tan26>=^ r . (9) 

m u ~ m s 

The charged squirk masses remain unchanged with mass 
m± = M v . 

Due to infracolor confinement, the squirks remain in a 
bound state. The dynamics of the squirk bound states 
will occupy much of the later discussion of the paper. 
We work in the approximation that cti C {m q ) < 1, so 
that the masses of the bound states are dominated by the 
constituent squirk masses. There is a hierarchy of mesons 
formed from orbital excitations of the bound squirks. We 
will be concerned mainly with just the J = L = (with 
necessarily S = for squirks) ground states, which we 
write as r)ij with ij labeling the constituent squirks. The 
charged bound states consist of 

V±o ~ (M±?4)>( u o«t) 



V±i ~ (u±u{), (umL) 



(10) 
(11) 

where the charged meson masses are roughly m^ ±j ~ 
m± + TOj . The neutral states consist of several distinct 
mass eigenstates 

7700 ~ (uquI) (12) 
V+- ~ (u-ul) . (13) 

?7u - (uml) (14) 



with masses that are again roughly m Vaa ~ 2mo, m ?7+ _ ~ 
2m±, and m mi ~ 2m\. Whether the heavier squirks 
hadronizc before weak decay, like the 6-quark of the SM, 
or decay before hadronization, like the i-quark of the SM, 
depends on the parameters of the model. As we will see, 
both regimes are relevant to the model. 



III. SQUIRK PRODUCTION 

Squirks can be produced through color and electroweak 
production. Conservation of infracolor implies squirks 
are always produced in pairs. At hadron colliders, squirks 
are generically produced with some kinetic energy for 
each squirk. As vividly demonstrated in [7], squirks will 
initially fly apart forming an infracolor string between 
them. The string will stretch without tearing until the 
energy in the infracolor string matches the total kinetic 
energy. The infracolor scale of interest to us is roughly 
Aqcd ^$ Ai C <C ^squirk, where the infracolor strings are 
microscopic, and thus the squirks shed their kinetic en- 
ergy and annihilate on timescales short compared with 
the timing systems of the collider detectors. 

Here we consider the several possible combinations of 
squirk pair production in terms of the squirk mass eigen- 
states. 



U±Uq + UQU 



The production of charged plus neutral squirks pro- 
ceeds through a s-channel W ± . This cross section is 
enhanced by the number of QCD colors and infracolors, 
while suppressed by the mixing angle associated with the 
U3 component of uq. The cross section at the Tevatron 
is shown in Fig. [2j incorporating the QCD color factor 
enhancement, but without the infracolor enhancement as 
well as without the neutral state mixing angle suppres- 
sion, to remain as general as possible at this point. 

Once squirks are produced, numerous effects must 
be considered to obtain a realistic estimate of the sig- 
nal rates. This includes understanding the spin-down 
process, /3-decay, and the various annihilation chan- 
nels. In addition, signal efficiencies are non-trivial when 
rn± — mo ^ Mw and thus the W decay is off-shell. This 
case, which is relevant for our model, implies the lepton 
Pt cuts tend to slice away more signal while also changing 
the dijet kinematics. A full estimation of our signal effi- 
ciency is not possible given the unknowns regarding the 
quirky spin-down dynamics. Instead, we have performed 
estimates of signal efficiency by constructing a "stand in" 
model, similar to Ref. [3j. This is discussed in Sec. VII 



B. UoUq and u\u\ 

The electroweak quantum numbers of Uo and U\ (T 3 = 
0, Q = 0) imply the couplings to 7 and Z exactly vanish, 
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FIG. 2. Leading order production cross section at Tevatron 
a(u±ul + uqu]^-). The three contours correspond to three 
choices of mo = 75, 80, 85 GeV from top to bottom. The cross 
section shown in this figure includes the QCD color factor 
and both electric charges, but does not include the infracolor 
(enhancement) factor, nor the mixing angle (suppression) as- 
sociated with the U$ component of uq. 



which has important implications. The first, and per- 
haps most important, is that uq and u± cannot be pair- 
produced at LEP II at tree-levelr] Second, the bound 

A J. ' — ' 

states formed from uqu and uiu[ cannot decay into pairs 
of photons or (potentially off-shell) Zs. 

The neutral squirks, uq and u±, can be pair-produced 
through QCD. The Tevatron production cross section 
of neutral squirk pairs is shown in Fig. [3] The lightest 
bound state of squirks must ultimately annihilate back 
into SM particles. There are several final state topologies 
that can result: gg and g ' g' (infragluon pair), illustrated 
in Fig. [3] as well as W + W~ (not shown). 

Although our model has negligible LEP II produc- 
tion of squirks, it is amusing to consider the size of the 
cross section for generalized squirks with (<3,T 3 ) arbi- 
trary. We find that cross section for colored squirks is 
roughly 0.6 pb for Q = 0, |T 3 | = 1/2. The cross section 
increases if Q ^ is taken. Since the total hadronic cross 
section at LEP II is measured to within ±0.5 pb to 95% 
CL |29j . even without considering the distinctive kine- 
matics (dijet invariant mass peak), we see that 80 GeV 
squirks transforming under the electroweak group are es- 
sentially ruled out. This is the motivation for our model 
employing an electroweak triplet that carries vanishing 
hypercharge. 



Charged squirk pair production, possible /3-decay, 
spin-down, and annihilation is a fascinating but intricate 



2 There is a set of one-loop diagrams involving virtual Ws and 
virtual u± exchange, but this is very small. 



FIG. 3. Leading order production cross sections at Tevatron 
for a(u±u\_) with (top curve; My — > m±) and a(uou ) or 

o(u\u\) (bottom curve; Mu — > n%o or m\). Again, like Fig. [2] 
the cross sections include the QCD color factor and both elec- 
tric charges, but do not include the infracolor factor. 



process. Charged squirk pairs can be produced through 
QCD, and so their production rate at hadron colliders 
is large. The cross section at the Tevatron is shown in 
Fig-i 

The signature of charged squirks depends crucially on 
the various competing timescales of spin-down, /3-decay, 
and annihilation. If both charged squirks were to rapidly 
/3-decay, u± — > W ± uo, this leads to a W + W^uqu' sig- 
nal, with UoUq annihilating into gg (or invisibly to g'g'). 
Whether this process is a feature or a constraint de- 
pends on several factors, particularly, the mass splitting 
m± — mo and the associated efficiency of detecting the 
W decay products. 

There is, however, another process that competes with 
/3-decay, namely the squirky spin-down and annihilation. 
As we will show in the next few sections, there is a region 
of parameter space where the spin-down and annihilation 
of u±u\_ occurs faster than /3-decay, leading to gg, g'g' , 
77, and W + W~ signatures. 



D. 



u±u 



The production of the charged plus the heavier neu- 
tral squirk proceeds through a s-channel W ± just like 
Sec. |III A| This cross section is enhanced by the number 
of QCD colors and infracolors, while suppressed by the 
mixing angle of the C/3 component of u\. While this 
process has a substantially smaller cross section than 
the one with ito, there are several interesting signatures. 
First, if the /3-decay process u x — > W ± u T occurs quickly, 
then this p rocess contains the same rich phenomenol- 
ogy of Sec. Ill C along with an additional W ± , which 
itself could be off-shell. If instead spin-down occurs 
rapidly, then a competition is set up between /3-decay 
with subsequent annihilation of 77^ , versus direct an- 
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FIG. 4. One of several diagrams illustrating squirky pair 
production through QCD of the lightest electrically-neutral 
squirks pp — > uou Q and annihilation «o^o ~* 99- Radiation 
(into color gluons, infracolor glueballs, etc.) is not shown. 



nihilation r] + i — > Wj. It is also possible that /3-decay 
of u± —> W ± Uo occurs, forcing ui — > W ± u T followed 
by u T — > W T uo. This gives a spectacular three- W plus 
2-jet signal. 

As we will see in the next several sections, the produc- 
tion of uqUq, u\u\, and u±u± are most relevant for the 
Wjj signal. 



given by 

r>( M± -> w±uo) = jV/G ^^ (m t' mo ' Wvy) (15) 

48-7r d 

where the function R is presented in its full glory in Ap- 
pendixlBl The kinematic function R scales approximately 
as | Amp, as expected for a 3-body decay. The mixing 
angle is defined by the transformation given in Eq.(d). 
The rate scales as the number of SM flavors and colors 
to decay into, Nt, which is 9 for m T < |Am| < m t + m^. 
The rate for u\ — > W ± u T is analogous, 

-> W^ u± ) = NfGWoR(^ ; rn ± , mw ) 

where the relative factor of 2 (compared with the charged 
squirk case) accounts for the two possible electric charge 
combinations. The inverse decay rate is roughly 

/o „ 2i \ /40 GcV\ 5 /0.717\ 2 

tf),u ± ~ 3 X 1CT 21 sec x — . 17 

V Ato J \ s e J 



IV. WEAK DECAY AND SPIN-DOWN 

We now consider the quirky dynamics that occurs after 
the squirk pair is produced. This includes /3-decay, spin- 
down, and the annihilation of the squirks into SM fields. 
The annihilation cross section itself is a sensitive decreas- 
ing function of the angular momentum of the squirk sys- 
tem. The angular momentum of the system is of order 
L 1 in the hard production, and then grows as the 
excited state losses energy through radiation. Since ev- 
ery radiated quantum is expected to change the angular 
momentum by AL ~ ±1, L grows roughly as the square 
root of the number of emitted quanta. As a result, the 
squirky bound state typically loses its excitation energy 
before annihilation [7j. On the other hand, /3-decay can 
occur during the energy loss process, as it is independent 
of angular momentum. Thus, in this section we discuss 
the competition between /3-decay and energy loss, and 
then in Sec. [V] we calculate the decay rates of squirko- 
nium states formed only after energy loss has occurred. 



A. /3-decay 

Weak decay of bound squirks or quirks can be readily 
approximated by considering the decay of an individual 
squirk in isolation. This is entirely analogous to heavy 
quark decay of (heavy) quarkonia in the SM. Charged 
and heavy neutral squirks can /3-decay through u± — > 
W ± uq and u\ — > W ± u T . When the mass splitting among 
these states Am^ = \m± — mj| is smaller than Mw, which 
is the main parameter region of interest to us, the rate is 



B. Spin-down 

The "spin-down" or energy loss process for squirks 
transforming under both QCD and infracolor is an in- 
teresting but difficult problem. In general, the energy 
loss can proceed through the radiation of a) regular QCD 
hadrons - mostly charged and neutral pions, or b) infra- 
color glueballs|^] Given the intrinsically non-perturbative 
nature of this emission, we will not assert which of these 
two distinct processes dominates, but instead simply 
make estimates assuming one or the other does, and dis- 
cuss the phenomenology. 

One reasonable approach to the mechanism of energy 
loss is to assume that at every squirky oscillation, the 
two QCD-hadronized squirks at the end of the infracolor 
string collide at semi-relativistic velocities, emitting a 
soft hadron, carrying energy roughly of order Aqcd or 
A ic [7J [5T] . The period of the squirk oscillations is of or- 
der n/Af c , where [i is the reduced mass of the two-squirk 
system. 

Consider first the case where only QCD dominates and 
the excited squirky system loses energy of order AE ~ 
Aqcd ~ 1 GeV after every oscillationFj The time to lose 
the total excitation energy Eij = Vs — (m, + rrij) for 



3 Photons may also play an important role in the case where 
quirks or squirks are charged and uncolored, but this will be 
subdominant to QCD radiation for our model. 

4 It is not unreasonable for the energy loss process to be domi- 
nated by QCD due to kinematics. This is because QCD contains 
pions, which are lighter than glueballs, unlike the situation with 
infracolor. 
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= 0,1, ±) is roughly 
fx Eij 



^loss j 



Af c Aqcd 

/ fJ- \ ( lGeV \ 2 f Ejj 
V50 GeV/ V A lc / ^100 GeV 
x (3 x KT 21 sec) . 



(18) 
1 GeV\ 



A 



QCD 



) 



For the case (i, j) = (±,0), the reduced mass /i = 48 — > 
53 GeV for m± 120 -> 160 GeV. Thus, we see that 
a typical excited squirkonium system settles down to its 
ground state on a timescale roughly comparable to j3- 
decay, Eq. (fiTl), when |Am| ~ 40 GeV. We stress that 
this mass difference is a rough approximation, given that 
we cannot determine the energy loss rate precisely. In 
this scenario, every squirky event will be accompanied 
by a large multiplicity of soft pions which will contribute 
to the underlying event. Such new contributions to un- 
derlying event physics were studied in JTTJ [5T] . 

Next, consider the case where the energy is lost also 
through infracolor glueball emission. This emission pro- 
ceeds in addition to the QCD emission described above, 
and so the timescale in Eq. ( p~8| ) is expected to be an 
upper bound. Infracolor glueballs generically have long 
lifetimes [7J, and thus will escape the detector as miss- 
ing energy. If infracolor glueball emission does dominate 
the spin-down process, then there is additional missing 
energy in every squirky event. This affects the px dis- 
tribution and particularly the T^t distributions of squirk 
signals, but we will not study this effect here (for a dif- 
ferent study on the effects of hidden radiation see [IE])- 



C. Survival Probability 

We can be more precise about our estimate for the frac- 
tion of squirk production that /3-decays before settling 
into the ground state r\. This calculation involves con- 
voluting the differential distribution daij/dEij with the 
survival probability, P(Eij) = 1 — exp[— t\ oss (Eij) /Uj] to 
obtain branching fractions: 



BR(u ± mJ ?/ ±0 ) 



1 



dE ±0 ^^P(E ±0 ) (19) 



BR(u ± ?4 -> W ± u ul) = 1 - BR(m ± 4 -» ?7 ± o) , (20) 
where i±o = tp,u±i an d 

BR(n + u + t r;+_) = — / ( 
o-+- J 



,IE-^-P(E + .) 



dE, 



BR(w +M+ t -> W + W~u ul) = 1 - BR(u +U+ t 



(21) 



(22) 

where f_| = tp }U± /2, where the factor of 1/2 accounts 

for /3-decay of either squirk. Lastly, for completeness, we 



emphasize that the lightest squirks cannot /3-decay, and 
thus lose energy to the ground state with unity probabil- 
ity, i.e., 



BR(uquJ) -> ?7oo) = 1 



(23) 



Here this branching ratio assumes the squirks lose energy 
before they annihilate, which is our working assumption 
for this paper [TT] . We use these results for our signal 
estimates in Sec. I VII 



V. ANNIHILATION 

Squirk pairs will eventually settle into a ground state, 
?7, which itself decays through constituent squirk annihi- 
lation analogous to quark annihilation leading to meson 
decay in QCD. 

In this section we compute the decay rates of the var- 
ious squirkonium states into SM fields and infragluons. 
The formalism for calculating the annihilation rates is 
similar to that of calculating the decay rate of quarko- 
nium in QCD |30j which has been adapted to quirkonium 
in [71 HI [TUJ [5D] . Moreover, squirkonium has a close rela- 
tive in supersymmetry, "stoponium" , whose decay rates 
have computed [3iTl33j . and we use to cross-check our 
own results. Under the assumption that the squirks set- 
tle into the ground state, ij (J PC = h ), the decay of 
7] — > ff is negligible, and will not be considered for the 
remainder of the paper. 

The squirkonium decay rate is 



r~au rel |^(0)| 5 



(24) 



where a is the annihilation cross section in question, v rc \ 
is the relative velocity among the quirks and | "0(0)1 is the 
squared wave function of the squirkonium bound state 
evaluated at zero squirk separation. The annihilation flux 
<7t> rc i is evaluated near threshold, where the two squirks 
are nearly at rest. Assuming the squirk wave function is 
Coulombic, the squared wave function is 



\m\ 2 = c v ^- 

TT 

where ap is the Bohr radius given by 



(25) 



a = [(C 2 (3)a s (a^) + C 2 CN lc )a lc {a^)) fi] 1 . (26) 

Note that there may be significant deviations from the 
Coulombic result as well as other QCD and infracolor 
effects. This can be partly parameterized by a coeffi- 
cient, c v , as discussed in Ref. [IB] . In practice, higher 
order QCD and infracolor effects can significantly mod- 
ify the Coulombic estimate, as can be found by scaling 
the QCD results of [51] to include infracolor as well as 
color. An additional assumption in our calculations is 
that we assume the Ni C dependence of the wavefunction 
associated with the initial bound state of a pair of squirks 
is to be treated the same way as the N c dependence of 
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QCD for quarkonia. In any case, we do not attempt to 
model the squirkonium potential beyond the Coulombic 
approximation, and thus take c n — 1 for our estimates of 
the absolute widths of r\. 

Now we consider the several possible annihilation chan- 
nels for the squirky mesons according to the constituent 
squirk states 

A. 7)00, 1711 



the infragluon emission. The infraglueball decay rate is 
generally extremely small, since it proceeds through a 
dimension-8 operator suppressed by squirk masses [7]. 
Infraglueballs are thus expected to escape the detector 
before they decay. In this case, the g'g' signal is itself 
invisible. There can also be a gg'g' process, when the 
squirks are produced in association with an additional 
initial- or final-state gluon (or quark) jet, leading to a 
mono jet signal. 



All of the neutral mesons can decay into QCD gluons 
and infracolor gluons, with decay rates 



1%. 

r(i7n g'g 



gg) 



3m? 
MNf c - 



\m\ 2 

-l)a 2 



2N ir 



mow 



(27) 
(28) 



for i = 0, 1. These rates are parton-level approximations 
without QCD or infracolor hadronization. In the limit 
Aqcd, Ai c <C too, the effects of hadronization on the de- 
cay widths is small [35]. Since Ui does not carry electric 
charge, the tree-level width of r\a — > 77 vanishes. There 
is also potentially a tree-level rate into W + W~ through 
a t/u-channel u±. This rate is suppressed by kinemat- 
ics (particularly for 7700 — > W + W~ , since m ~ M\y) 
as well as the weak couplings associated with this chan- 
nel. While potentially interesting for r/n, which can be 
copiously produced at LHC, it is subleading compared 
with the above glue and infraglue rates, and we do not 
consider it further. 

The gg channel can dominate, allowing uqUq produc- 
tion at a hadron collider to lead to a dijet bump. For 
the Tevatron, the production rate is far smaller than the 
QCD background, leading to neither an observable sig- 
nal nor a constraint [27]. Interestingly, UA2 has com- 
paratively strong bounds on dijet resonances. The pro- 
duction cross section cr(u ?4) a * UA2 (y/s = 630 GeV) 
with rn — 80 GeV is less than 1 pb per infracolor. As- 
suming uoUq quickly spins down and annihilates into a 
dijet resonance of about 160 GeV, the limit we should 
compare to is the W /Z' limit at the same mass. This 
is a (dijet) < 90 pb for M w , / z , — 160 GeV, and so is 
completely safe. 

The next comparable channel is annihilation into in- 
fragluons: g'g'. The signature of the g'g' topology de- 
pends on the decay rate of the infraglueballs formed from 



5 There is also the possibility that quirks or squirks could form 
"hybrid mesons" in which the quirk pair is in an infracolor non- 
singlet state 25 , leading to a set of qualitatively different decay 
processes. These include the single emission of glue or infraglue, 
combined with another SM gauge boson. The dynamics of in- 
fracolor non-singlets is very difficult to estimate without a more 
complete picture of the quirky hadronization dynamics, and so 
we do not consider it further. 



B. rj± 

The r)±o represents an electrically charged squirk plus 
neutral squirk system that has settled into its ground 
state before the u± has had a chance to /3-decay. The 
dominant annihilation channel is r]± — > W 7, with rate 

m 



±0 



-Kaa w N c N ic sl 
(toq + TO ± ) 2 



a/sl 



f(m ,m±,m w ) \4>(ty\ 2 
(29) 

and the kinematic function 



where a w — <^/o w , 
f(mo,m±,mw) has been relegated to the Appendix |B| 
Whether the timescale for this rate is slower or faster 
than /3-decay depends on the parameters in the model. 
Small mass splittings tend to suppress /3-decay, whereas 
small Ai C , which leads to a larger Bohr radius, tends 
to suppress the annihilation rate. There are potentially 
additional decay channels, including WZ and Wh, but 
these are much more suppressed by phase space, and so 
we will not calculate their widths. 



V+- 



The neutral meson formed from the two charged 
squirks can also decay into QCD gluons and infracolor 
gluons, with decay rates identical to Eqs. ( 27|28 ) sub- 
stituting m, — > m±. Since the constituent squirks are 
electrically charged, r}-\ can also decay into into pho- 
tons with a decay rate 



77)= 6 ^r L W0)| 2 



(30) 



This is an interesting annihilation channel, and poten- 
tially a constraint on our model, given existing bounds 
on 77 resonances from collider data. We discuss the size 
of this signal in the next section. 

Finally, there is another uniquely squirky annihilation 
mode, namely rj^ — » uo^o, through the quartic inter- 
action, Eq. ([2]), proportional to A4. The rate for this 
annihilation is given by 



u u\) 



9TTNl(a iS 2 e + a U4 c 2 9 ) 2 ) 



2m| 



#|V(0M&) 
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??±0 




/3-decay 



W ± U ul 




/3-decay 



W ± ~f W ± u ul 



FIG. 5. The "tree of life" for the charged + neutral squirk 
pair. While only u±u is shown, the same tree of life also 
applies to uqu^. 



u±u± 



loss 




/3-decay 



V+- 




/3-decay 



99 77 9 9 



UqUq 



u ± u' 



FIG. 6. The "tree of life" for the charged squirk pair. Charged 
squirk pair tree of life. Everywhere u±u\ j is written is under- 
stood to mean u±u\ ) + uqu^. 



where a^ui = ^4^4/(4^). This so-called "double 
wonga-wonga" process is unique in that the annihilation 
is from heavy squirks to light squirks without any hard 
SM emission. In practice, the light squirks will subse- 
quently spin-down, analogous to the spin-down phase fol- 



lowing ordinary collider production of squirks, and then 
the light squirks bind up into an 7700 and annihilate as 
discussed in Sec. IV Al 



VI. BRANCHING FRACTIONS AND SIGNAL 
ESTIMATES 

With the various decay rates in hand we can now com- 
pute the total cross sections for several final states of 
interest. We will again separate the discussion according 
to the different initial states. 

A. u±ul + uqu^ 

The detailed branching fraction "tree of life" of the 
quirk pair u±Uq + uqu^ is shown in Fig. [5] The first 
branch splits between /3-decay of u± (right branch) versus 
radiative energy loss, leading to the ground state meson 
77-1-0 (left branch). The relative probabilities of these two 
branches are set by Eqs. (19 1 and (20). 



Following the left branch, in which the ground state 
77-1-0 is reached, the tree of life then branches further, as 
the meson 77-1-0 can decay into several different channels. 
Considering all possible /3-decay branchings of u±Uq into 
W ± u (j u\ )1 the relevant branching fraction here is 



BR(t7 



±0 



W ± u ul) = 



W-y 



(32) 



where the Tp is given by Eq. ( 15 1 and T(Wj) is given by 



Eq. (29) 



Assuming one of the /3-decay branches is taken, the 
tree of life results in the squirk pair uqUq, which after 
energy loss, ends in the ground state 7700- The branching 
fraction of 7700 into visible SM particles is given by 



BR(t7o 



99) 



r(r?oo -> 99) 



r(??oo -> 99) + r ( ? 7i 



00 



(33) 



where the decay rates are given in Sec. |VB| 

Combining these branching fractions, we obtain the 
cross section of u±u\ + u Q uL into the Wjj final state, 



o±o{Wjj) = a(pp -> u±ul) BR(u±ul -> W ± u Q ul) + BR(u±uJ -> ?7±o)BR(77 ± o -> W ± u u\) BR(t7 00 -> gg) (34) 



where the two terms in the square bracket represent the two paths to our final state, determined whether /3-decay or 
energy loss occurred first. The cross section into the Wj final state is 



a ±0 (Wj) = <r(pp -> u ± uj) BR(tt±uJ, -> 77 ±0 )BR(77 ±0 -> W ± -f) 



(35) 



This final state has a bound of about aw-y 8-13 pb to 
95% CL, using constraints from DO [35] and CDF [37]. 
The cross sections given for a(pp u±Uq) are under- 



stood to sum both u±u\ and uqu^- contributions. 
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u±u± 



The tree of life for the u±u± is somewhat more com- 



plicated and is the result of combining Figs. |6|5| The 
branching begins between energy loss and /3-decay (which 
is now twice as fast to account for the decay of either 
squirk) at the top of Fig. [6] If energy loss is fast, the sys- 
tem will typically get to the ground state meson ry + _ (left 
branch). Here there are several interesting decay modes, 
including gg, g'g' (invisible), 77^] and uqUq. The latter 
can play an important role in cases where there is tension 
with the 77 or pre-tagged top constraints. In particular, 
if A4 ~ 2, the mode Eq. (31 ) can dominate, effectively di- 



We note however that diphoton and pre-tagged top are 
safe in large regions of our parameter space even in the 
absence of the A4 coupling. 

The branching fraction of 77-1 to /3-decay is 



BR(r ?+ _ -> W T u±ul) = 



2T f 



minishing the branching fraction into dangerous modes. 



2r /3 + T gg + r sV + r „ „t + T yl 

(36) 

where u±Uq is understood to mean both u±u and uou^. 

We pay close attention to the WWjj final state, which 
is the result of the /3-decay of both squirks and the subse- 
quent annihilation of 7700 into dijets. Combining branch- 
ing fractions, we obtain the cross section for the "two- 
armed lobster" shown in Fig. [T] into the WWjj final 
state, 



o+-(WWjj) = a{pp —> u±u< ± ) BR(u±4 W T u ± ul) + BR(w ± M f ± 77 + _)BR(r7 + _ ->■ W T u±u\ 



BR(u±itQ 



W ± u ul) 



BR(u±t4 ->• r7±o)BR(r/ ±0 -»• W ± u ul) BR^oo -»• gg) (37) 



where the two terms in each square bracket represent the four distinct paths to our final state, determined whether 
/3-decay or energy loss occurred first in either of Fig. [|or[6j Again, u±u\ is understood to mean both u±u\ ) and uqu^- 
in the branching ratios, as appropriate. The cross sections of u±u± directly into various resonance final states are 



= a (PP -> w±«±) BR(it ± i4 -> 77 + _)BR(?7 + _ -> ij) 



(38) 



where ij can be any of gg, g'g' (invisible), 77, and wotij. 

The WWjj final state is important for two reasons. 
First, the WWjj signature can "leak" into the Wjj 
search when the decay products of one W are either lost 
or do not pass the CDF analysis cuts for their exclusive 
analysis. In particular, if the mass splitting m± — mo is 
not too large the jets from a hadronically decaying W 
will be softer, and can frequently fall below analysis cuts 
on jet Et- The CDF collaboration analysis for Wjj has 
a jet energy requirement of Et > 30 GeV, suggesting it 
is much more likely for the jets from a hadronically de- 
caying off-shell W to not pass their Et cuts. This means 
that charged squirk pair production provides an addi- 
tional production source of the Wjj signal after detector 
cuts. As the mass splitting is increased, a decreasing 
fraction of jets from a hadronically decaying W would 
fail their Et cut, leading to a smaller contribution to 
the Wjj signal. 

Second, consider the dileptonic process, where both 
W's decay to e or ^. The signature of this process, 
11 jj plus missing energy, is nearly identical to top 
quark production and dileptonic decay, tt — > W + W~bb, 
when 6-tagging is not required of the jets. This was an- 



Note that the final states W^W^ and ZZ may also be of inter- 
est, particularly at LHC, but are not calculated here. 



alyzed by CDF, where their measurements of the fully 
leptonic "pre-tagged" top cross section implies an upper 
limit on this signature of about 2 pb at 95% CL 38 . 

This is likely a reasonable limit for on-shell VF-pair 
production in association with 2 jets. However, off-shell 
Ws can lead to suppression of this signal due to the 
smaller fraction of events with sufficient energy to pass 
the detector cuts. We make some rough estimates of the 
efficiencies of both the Wjj and pre-tagged top analysis 
in the next section. 

Another final state which can potentially constrain the 

parameter space is the decay of the 77^ meson into 

77. Searches for diphoton resonances were performed 
by CDF [3S], DO 0D], and CMS gT] in the context of RS 
graviton searches. The constraint on the cross section for 
a diphoton resonance, in the mass window of interest to 
us, 2m± ~ 240-290 GeV, is of order 10-40 fbQ 



7 A more precise number is not straightforward to extract because 
CDF 39 has a diphoton excess at an invariant mass of 200 GeV, 
leading to a weak bound of 100 fb (as opposed to the expected 
30 fb), while DO I40| did not present an exclusion plot for 77 
resonance cross section independent of the ee resonance cross 
section. 



VII. BENCHMARKS AND EFFICIENCIES 

We present two benchmarks with the parameters, 
masses, cross sections, branching fractions in Table [IT] 
The two benchmarks represent two qualitatively differ- 
ent regimes: Benchmark 1 has rapid spin-down, and 
then squirky /3-decay or annihilation, whereas Bench- 
mark 2 has rapid /3-decay, and then annihilation. Each 
Benchmark results in qualitatively distinct signal kine- 
matics, as we will explain. For each Benchmark, we also 
present our estimates of the relevant limits on certain 
parameters and signal rates. The input parameters are 
N ic ,My,Mg,5,X4. We took N ic = 3,4 infracolors for 
Benchmark 1,2, leading to a signal rate into dijets that 
easily satisfies the UA2 bound. We took My , Ms, 5 such 
that the masses worked out to m±,mo — 120,80 GeV 
and m±,m = 145,75 GeV. This allows us to illustrate 
the qualitative differences between Am = 40 GeV ver- 
sus Am = 70 GeV. We chose Ms to be slightly heav- 
ier than My, such that the contributions to the isospin- 
violating electroweak correction, AT, essentially vanish, 
as described in Appendix [X] Finally, we took A4 to be 

of order one, which ensures the cross section (77) is 

less than 10 fb. Smaller values of A4 imply larger rates, 
of order tens of fb, into diphotons. 

The cross sections cr(uoi4), ct(u±Uq + uqu^-), and 
a(u±u±) can be read off from Figs . [2I3I Weak production 
is suppressed by couplings, while the colored production 
of u{u±v} ± ) is enhanced by couplings but further sup- 
pressed by kinematics, leading to rates at the Tevatron 
that are roughly comparable. 

Each model has A^ c <C m squ i r k, such that the infracolor 
coupling, a ic (m g ), is perturbative when evaluated at the 
mass of the squirk. The choice A^ c ~ few GeV implies 
infraglueballs decay well outside the detector, but de- 
cay fast enough to not cause cosmological conundrums. 
Hence, the g'g' final state leads to no hard SM parti- 
cles, and possibly missing energy (depending on the de- 
cay kinematics). 

The cross sections a±o(Wjj) and a + -(WWjj) corre- 
spond to Eqs. (34) and (37) respectively. These cross 



sections form the starting point for obtaining the Wjj 
signal, as well as several additional signals for which the 
Tevatron has already placed constraints. 

We have then estimated the efficiency to detection 
for three signals: the efficiency for a±o(Wjj) to pass 
the CDF Wjj analysis cuts [T], the efficiency for 

cr_| (WWjj) to pass the CDF Wjj analysis cuts (where 

one Ws decay products are missed or not energetic 
enough to pass the CDF cuts), and finally, the efficiency 
for a + _ (WWjj) to pass the CDF ti pre-tag analysis cuts 
[3S], which we call a(l + l-^jj) x cff. 

Our estimates are based on a "stand-in" model for 
squirk production and decay, modeled as intermediate 

resonances with masses the same as , r\±§ and 7700; 

allowing for /3-decay into r\±§ which in turn can /3-decay 
into 7700- We have implemented both of these models in 
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0.49 
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^ 90 


cr±a(Wjj) 
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a (Wjj) x eff 


1.3-2.0 
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U.UUD 


U.UU4 


^ U.U1-U.U4 


o ±0 (W~i) 


1.1 


0.2 


< 8-14 


AT 


0.02 


0.01 


-0.05 -> 0.2 


0"LHC7(^O^o) 


480 


430 




0"LHC7(M±M±) 


200 


130 





TABLE II. Benchmark models in parameter space. All masses 
in GeV, all cross sections are in pb for Tevatron (except 
"LHC7" for yfs = 7 TeV LHC and "UA2" for yfs = 630 GeV 
CERN SppS). The cross sections are discussed in Sec. |III| 
the branching fractions into various final states discussed in 
Sec. |VI| the efficiencies are discussed in Sec. |VII[ and AT 
calculation is done in Sec. [X] 



MadGraph [42] with various mass splittings. The model- 
ing of missing energy cannot be reliably done, given the 
squirky spin-down process that can emit infraglueballs 
which escape the detector as missing energy. 

Nevertheless, we can implement the various lepton, jet, 
missing energy and transverse mass cuts on the "stand 
in" model to obtain "stand in" efficiencies which should 
give us a reasonable idea of where the squirk model stands 
with respect to the various signals after cuts. Both 
benchmarks lead to some excess in the dileptonic pre- 
tag ti sample (a(£ + £~( )jj) plus missing energy), but 
within the 95% CL limit from CDF. Both benchmarks 
also lead to between 1-2 pb of Wjj signal, where the 
range corresponds to including (not including) the miss- 
ing energy and transverse mass cuts for the lower (upper) 
end of the range shown. Also shown in Table [TT] is the 
ratio WWjj/Wjjtotai which corresponds to the fraction 
of Wjj signal after efficiencies that arise from the two- 
armed lobster versus the one-armed lobster. We conclude 
that the squirk model we have presented is capable of 
generating a Wjj signal consistent with the CDF excess, 
so long as the efficiency of detection of squirks is compa- 
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rable to our "stand in" model for event-level simulation. 

While we have presented estimates for a variety of 
observables above, we emphasize that several quantities 
should be taken as rough estimates due to the various un- 
certainties involving quirky dynamics. For example, the 
uncertainty on the timescale of energy loss can strongly 
affect the branching fractions between spin-down and /3- 
decay as shown in Fig. [5] and [6} In addition, squirks that 
rapidly lose energy often reach the r\ ground states, which 
have definite masses, and so leads to features in kinematic 
distributions. Even when energy loss is rapid, the choice 
between branchings involving /3-decay (such as Wjj and 
WWjj) versus annihilation channels (such as Wj and 
77) is sensitive to the uncertainty in the wavefunction at 
the origin for quirkonium states. 

VIII. IMPLICATIONS AND DISCUSSION 

We have demonstrated that the Wjj signal observed 
by CDF can be obtained from a model of squirks without 
violating existing collider bounds. There are two quali- 
tative regimes where the signal arises: 

1. Rapid spin-down, then squirky /3-decay or annihi- 
lation. In this scenario, the excess energy squirks 
carry, — 2m, is released quickly, allowing the 
squirk pair production to settle into squirkonium. 
This could happen because /3-decay is suppressed, 
or could happen if the energy loss is rapid, or both. 

A competition is set up between /3-decay of the 
squirkonium constituent squirks versus direct an- 
nihilation of the squirks. In this regime, several in- 
teresting squirkonium decay modes with invariant 
mass resonances are potential confirmation signals, 
including gg, Wj, 77, g'g' (invisible), and possi- 
bly W ± W^, W ± Z, ZZ, and signals with a Higgs 
boson. 

2. Rapid /3-decay, then spin-down, then annihilation. 
In this regime, the only squirkonium state that 
is reached is 7700, yielding the famous jet-jet res- 
onance. This could happen if /3-decay is compara- 
tively fast, or could happen if energy loss is com- 
paratively slow, or both. 

In this regime, every heavy squirk /3-decay yields 
a (possibly off-shell) W, giving many signals with 
multi-W^'s plus a jet-jet resonance. 

The two regimes share several interesting signals in- 
cluding the Wjj final state. Interestingly, neither regime 
contains an associated "fjj or Zjj signal, which is char- 
acteristic of this model. We emphasize that while squirk 
production and squirk /3-dccay rates can be calculated 
perturbatively, the "spin-down" or energy loss phase, as 
well as the squirkonium decay rates have substantial un- 
certainties. The relative branching fractions of squirko- 
nium can in some cases be determined, since the decay 



rate dependence on the wave function at the origin drops 
out. But determining which of energy loss or /3-decay is 
more rapid can merely be estimated. 

The kinematics of the two regimes to Wjj are quali- 
tatively distinct. In particular, if the second term dom- 
inates one would expect an invariant mass peak for the 
whole Wjj system, and an edge in the transverse mass 
distribution, while if the second dominates such features 
would be either absent, or less pronounced. CDF have 
presented interesting distributions [2J, but we leave the 
investigation of such kinematic features for future work. 

We have calculated the LHC production cross sec- 
tions of the squirks in the parameter ranges relevant to 
this model. We find that the electroweak production of 
charged plus neutral squirks is relatively small, typically 
a few pb. The colored production of squirks is, not sur- 
prising, quite large - of order hundreds of pb! The LHC 
production of u±u± seems particularly important, since 
depending on which regime we are in, it can yield sig- 
nal rates into the annihilation modes [regime 1] or the 
multi-lF production [regime 2], 

There are several associated signals of our model that 
we have not discussed. Squirks have nontrivial interac- 
tions with the Higgs boson through the dimcnsion-4 op- 
erator, Eq. ([3]). Loops of squirks modify the effective 
hgg coupling as well as /177. Since Uoi u i receives its 
mass from both electroweak preserving, Eq. (JlJ and elec- 
troweak breaking Eq. ((3]), the contribution will scale as 
S 2 /(6 2 + M 2 ). 

In this paper we have focused mainly on a quirky ex- 
planation for the CDF Wjj excess, but in the process 
we have illustrated many interesting signals of quirk or 
squirk production at colliders. With the Tevatron hints 
of new physics, combined with the wonderful prospects 
at LHC, we expect an exciting time in quirky physics. 

Note added: Just before this paper was completed, 
DO reported an analysis of their Wjj data, finding "no 
evidence for anomalous resonant dijet production" [43] . 
Taken at face value, their analysis sets an upper bound of 
1.9 pb at 95% CL for resonant dijet production near the 
invariant mass window of the CDF excess, and thus does 
not rule out a new physics explanation of the CDF excess 
with a cross section less than this value. Given that the 
DO analysis was done with lower luminosity (4.3 fb _1 [13] 
versus 7.3 fb _1 [2]), without an inclusive (2 or more jets) 
analysis (unlike [2]), we remain skeptical of drawing neg- 
ative conclusions until the joint task force |44j completes 
their analysis. 



Appendix A: Electroweak Precision Corrections 

The introduction of the Higgs operator at dimension-4, 
Eq. (J3J) , splits the masses of the fields within the triplet 
V. This isospin violation leads to modifications to elec- 
troweak precision data. The correction to 5, which char- 
acterizes Bp O Wp mixing induced by electroweak sym- 
metry breaking, vanishes here since the squirks transform 
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mV [GeV] 



FIG. 7. Contours of AT per infracolor, but including QCD 
color, in the (mu,ms) parameter space. From left to right the 
contours are -0.1, -0.05, -0.025, +0.025, 0.05, 0.1, 0.25. The 
"funnel" region has AT ~ due to a cancellation between 
the loops of the light and heavy neutral scalars. 



in electroweak representations with zero hypercharge. 
The correction to T arises from the mass differ- 



ence between the charged and neutral components of 
the isotriplets. This has been calculated for a general 
scalar multiplet with arbitrary isospin and hypercharge 
in Ref. [351 . Applying their results to our case, we obtain 



AT = 



N n Ni, 



[cjf(m ,m ± )-slf(m 1 ,m ± )] (Al) 



in terms of the self-energy loop function 

f ( \ 2,2 2m2 afn 2 b 

j{m a ,m b ) =m a + m b 



log- 



(A2) 



As expected for a renormalizable theory, the correction 
to T is finite, and vanishes in the various limits: too — > 
TOi (the operator vanishes); my ^> ms (decouple the V 
scalar); and ms 3> my (decouple the S scalar). 

Since the model contains negligible additional con- 
tributions to AS, the quantity AT can be bounded 
from fits to S and T, e.g. 46J, where one can allow 
-0.05 < AT < 0.2 and remain within the 95% CL limits 
for m/i = 115 GeV. Wide ranges of parameters allow 
for sizeable splitting between the charged and neutra l 
squirks (allowing weak decay to proceed, c.f. Sec IV A I, 
while having T easily within the electroweak bounds. 



Appendix B: Collected Formulae 



Here we collect some analytic formulae used in earlier parts of the paper. The function / used in Eq. (35) for the 
weak annihilation of r)±o into W7 is 

[(to + to ± ) 2 - m 2 w \ [m w (-3ml + 4m m± + 3m|) + to (to^ - m 2 ± )(3m + m±)] 

f{m ,m±,m w ) = ± ^ ^— — \ J - (Bl) 

2m Q m^.m w (m Q -I- m±) 

y rate of u± — > u + W* 

9m^y (m4 — TOq) + 2771^, (toq — to 2 ± ) 3 + 6to^(toq — m 2 ± ) 



The function R used in the formula for the beta decay rate of u± — > u + W* in Eq. ( 15 ) is 

')m'.. in)' — in':\ — ''hit, I hi- — m~ 
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where 



and 



m + m ± + m w - 2m w (m + to ± JJ log [ 



A 2 = \J (to + m± — mw)(m — m± + mw)(—mo + m± + tow)(to + m± + mw) 

A 2 = \J (—too + m± + m\y)(mo + m± — mw)(m a — m± + mw)(m + m± + mw) 
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